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Abstract 

. 

C/3 ! We determine the solution of the fractional spatial diffusion equation in n-dimensional 

£T) ' Euclidean space for a "free" particle by computing the corresponding propagator. We 

employ both the Hamiltonian and Lagrangian approaches which produce exact results for 
the case of jumps governed by symmetric stable Levy nights. 
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1 Introduction 



In nature there is a broad variety of systems in which the correlations in space or time give rise to 
anomalous transport whose probability density function (pdf) is non- Gaussian and the squared 
displacement is non-linear in time or diverges. A typical example is an asymmetric particle 
plume which spreads at a rate inconsistent with the classical model. We will concetrate on the 
power-law pattern EXl ~ K a t a in the super-diffusion case 1 < a < 2, which is manifested in 
a diverse number of systems £Q . 

In a continuous time random walk framework j2] if the jump length and the waiting time 
are considered to be independent random variables the joint pdf ip(x,t) can be written as a 
product of the marginal pdf 's for each random variable. It can be shown that if the waiting 
time satisfies a Poisson distribution while the jump length follows a symmetric stable Levy 
distribution (see Appendix A for the definition) then a Markov process (the first moment of 
the waiting time is finite) is generated and satisfies the fractional spatial diffusion equation jH] 

= K°T%P(x,t). (1) 

In (0) the real parameter a ranges into the interval (1, 2], the generalised diffusion constant has 
dimensions [K a ] = ^L- with [L], [T] being the dimensions of length and time, and D™ is the 
Weyl operator as defined in the Appendix A. The solution of (JIJ with the sharp initial condition 
lim 4 ^o+ P(x,i) = 5(x) can be obtained analytically in terms of the H-function jl] (the general 
definition of Fox functions is given in Appendix B) 



P(x, t; a) = — 7-rHH 
a\x\ 



\x\ 



(LI) 
(1,1), (1,|) 



(2) 



Its structure depends exclusively on the underlying geometry of the space. Taking in (J2J) the 
limQ.^2- P(x, t; a) the classical Gaussian solution is recovered by standard theorems of the Fox 
functions. 

The present work aims at the construction of the quantum mechanical path integral over 
Levy flights. The outline of the paper is as follows. 

In Section 2 the matrix elements of the time evolution operator in real time are calculated 
by inserting complete sets of eigenstates of momentum and position operators. Using the 
property that plane waves remain eigenstates of the Hermitian operator (— V 2 ) * and performing 
the integrations over the intermediate positions we end up with an analytic expression for 
the propagator of the "free" particle in Euclidean space in one- and three-dimensions. The 
propagators are then written in closed form in terms of the Fox functions. We also derive the 
asymptotic behaviour of the transition probability density in one dimension for large and short 
distances. In two- and n-dimensions (n > 3) the propagators are written in integral form using 
the Bessel functions of the first kind whose order is half an odd integer. 

In Section 3 we write the propagator in a path integral representation in which all paths 
contribute in the quantum mechanical evolution but they are weighted with the complex weight 
eri M , with Sm representing the action of the particle. Descritizing and decomposing the paths 
into classical and quantum fluctuating trajectories and using the binomial expansion we derive 
the asymptotic behaviour of the propagator in the long time limit at fixed position. Our result 
is then confirmed by employing the steepest descent method to the result extracted from the 
Hamiltonian approach. 
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In Appendix we give all the mathematical background needed for understanding some no- 
tions of the probability theory as well as integral relations which help to manipulate the cal- 
culations of the text. We also prove a proposition that justifies the use of plane waves as 
eigenfunctions for the operator (— V 2 ) 2. 



2 The Hamiltonian approach 

We consider the time independent, one-dimensional Hamiltonian operator 

H(P;a)=(j^y , ae(l,2) (3) 

and investigate the matrix elements of the time evolution operator , the so-called propa- 
gator 

P(x in ,Xf,j3;a) =< x f \e~^ H \x in > . (4) 

In (J3J \xi n > and < Xf \ stand for the initial and final states of the particle which are eigenstates 
of the position operator x 

•^\^in ^* •Ein\%in ^ Xj <C X j\. (5) 

The position and momentum eigenstates form a complete and orthonormal set of states 

J dx\x >< x\ = I = J dp\p >< p\ 
< x\y >= S(x — y); < p\p >= 5{p — p') (6) 



with inner product 



< x\p >= — -e^ p ' 3 (7) 



where the angle between p and x is zero or ir. We insert in (J3J) N complete sets of momentum 
eigenstates and N-l complete sets of position eigenstates. The time lapse between the initial 
and final states has been discretized according to (3 = Ne with e > 0. We then obtain 



N 



k=l 

N 



P(x ,x N ,p;a) = <x N \e ^ H \x >= ] \ dp k dx k ^ x < x k \e * H(p) \p k >< Pk\x k -i > 



1 r r 

= (^pliy J dpkdx,^ er^"-^-i H ^ (8) 

with the identifications \x in >= \x > and |x/ >= |x;v >• In (JHJ) we have used the property 
that plane waves are eigenfunctions of the Hamiltonian with eigenvalues \p\ a (see Appendix A) 
therefore 

e-i A ®\p k >=e-t H ^\p k > . (9) 
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Performing the integrations over x^s we obtain 



P(x ,x N ,(3; a) 
Using the identity 



n 



N-l 



2nfi 



N 



roo N N 

/ [ (dp k eiw^-^-O-f^tffc) Y[6(pi 

J -°° fc=l 1=2 



Pi-i). 



(10) 



5(p k - Pk-Je-t^-^Sfa-! - Pk - 2 )d Pk -i = 5( Pk - p k _ 2 )e~i H ^ (11) 

• , pi we arrive at the transition prob- 



and successively performing the integrations over P n-i 
ability density 



P(x N - x ,P; a) 



1 



2ixh 



j- PN (x N -x )-^H(p N )^p N _ J_ l°° e ip N (xN-xo)-cP(\p N \ 2 )% 



2lX 



~dp N . (12) 



(2m) TT 



will play the role 



In (I12|) the momentum has been rescalled and the parameter c 
of the diffusion constant. We recognize expression (J12J) to be the Fourier transform of the 
characteristic function of the stable law S„(0, — — n-, 0). 

QV ' (2m) 1 ' y 

One can prove that ()12j) is the solution of the fractional spatial diffusion equation 



dP(x,t) 
dt 



c(-V 2 ) 2 P(x,t). 



(13) 



To prove this we set for simplicity x = 0, x^ = x, P = t and partial integrate (JT^J) with respect 
to time in which case it reveals that 



dP(x,t) 
dt 



2tt 



roc a 

c e ipx \p\ a e- c(3{lpl ^dp= -cf- 1 



\p\ a e 



-c(3(\p\ 2 )? 



c -V 2 2 P(x,t).(U) 



In deriving the second equality in (fUjl we used the result ()A.7|) of Appendix A. 

One can proceed by writing the /3-dependent integrand of (|12[) in terms of the H-Fox function 

as 



ft I 2m 



1 



1,0 



— -"0,1 

a 



f py\PN 



2m 



a 



(15) 



thus having 

I(\x N - x \;p,a) 



2nh 



a\XN — Xq\ 



-H 1 ' 1 

n 2,2 



FAT — #0 



i,i 



(16) 



The propagator is then given by 

i 

P(x N - x ,P; a) 



2-Kh 



I(\x N - x \;P,a) 



1.1 
2.2 



f 2m (fi\ a 



- x N - x 



(l,y (1,1) 
(1,1) (1,1) 



(17) 



a result that was already discussed in [2] by solving (J1J) in Fourier-Laplace space. 

The asymptotic behaviour of the propagator for \x^ — x \ — > oo or \x^ — x \ — > at fixed 
time, can be determined by partial integrating (fT2*j) and rescaling the momentum leading to the 
expression 



Bca r°° ^ , ^ — £ — r V 



- x , (3; a) = — / sin £e V i^-ol ) ^. (18) 



7T|Xat — Xo| a+1 

• Large x expansion 

The propagator with the help of (fTHj) is 

P(x N -x ,P;a) « — — T T(a) sin(— ), a<2 (19) 

which exhibits the known power-law tail of the Levy distribution. Due to this property, 
the mean squared displacement diverges EXf — > oo. 

• Short x expansion 

Taylor expanding the cosine and changing the integration variable to w — /3c\p\ a in f)12j) 
we have 

P(x N - x , # a) « — — ^T(-). (20) 
7ra(/5c)^ « 

The generalization of the propagator in n-dimensions can be achieved by considering the 
Hamiltonian operator 

Hip; O (21) 
In this case the transition probability density is 

roc roc a 

P(\x N -x lP;a) = -—- •••/ dV e nPN<*»-*°H n<P») (22) 

' J 2m N ) anc ^ re P ea t e d indices are 
summed over. In deriving (|22|) we used the proposition of the Appendix A. 

It is obvious that expression depends only on the magnitude but not on the direction 
of x*n — x . We distinguish the following cases 

1. Two- dimensions 

We adopt plane polar coordinates, choose the first axis of the pV-space in the direction 
of xn — Xq, and find 

/*oo q p2n 

P(\x N -x Q \,p-,a) = — — ^ j o dp N p N e-K H ^J o dden^-^ cosd . (23) 

But J 27T d9 e lccosd = 2nJ (c), where J is the standard Bessel function defined in Appendix 
C. Thus 

P(\x N -x \,f3;a) = -^ f° dp N p N e ~^ H ^ J (^f\x N - x \). (24) 
2nTi Jo n 



5 



2. Three- dimensions 

Adopting spherical coordinates we find 

P(\x N -x \,/3;a) = 



dpN p 2 N e-^ H{PN) / (dcosfl) e iPN\s N -M^e 

(ziXTl) 6 JO J-l 



2 77 



d 



27T 2 /l |xtv — Xo| d|j?Ar — X | Jo 

1 1 d 

2% 2 h \x N — x \ d\x N — x \ 



ch\xn — %o\ 



1,1 

2,2 



'2m ( fi 



h \(3 



\x N - x \ 



(1.5) 
(1,1) 



(25) 



3. n- dimensions, n > 3 

Let da denotes the uniform probability measure on the unit sphere S*™" 1 = {x G i? n : 
||x|| 2 = 1}. Then dV = f^j^dp^da where the constant factor represents the area of the 
sphere in n-dimensions. The squares have a Dirichlet Di(^, ■••,!) joint distribution 
so each aj is distributed as the square root of a Be(^, 2 random variable. Thus 
the integral becomes 



P{\x N - x \,P; a) 



27T2 r°° 



(27rh) n r(f ) Jo Jsn-i 

1 27rt 

(27r^r(i)r(2=i)-/o 



d^da pSr 1 cos( ^- ( ^~ fo) )e-f^) 
dpv p A T 1 e-« H ^)F(p Ar , |^ - foUX26) 



where 



^(Piv, Wn -x \,n) = / 

Jo 



1 , /Pjv|a?jv - «o| 

aw cos( 

o n 



y/u)u 2(1 — u) V 



71 — 1. 



2fi 



n-2 
2 



2 \Pn\xn - x \ j 
Combining (}26|) and (|27|1 we finally have 



Jn^x{-p N \x N -So|). 



P(|fiv - ^o|,/3; a) 



1 



1 



(27T^) 2 ^|xat — Xq I 2 



7[G(\xn - x \,n) 



(27) 



(28) 



where G(|xat - x \,n) = J °° dp N pfae~ri H{pN J^- y {hp N \x N - x \). 



lr The pdf of the Dirichlet distribution with parameters (a±, • • • , a„+i) is 

ti \ r ( a i "I l-On+l) an -i I 1 

/(ii,'-,x n ;ai, ■■■|tt„ + i) = ^ ^ ^- -ij ■•■2C n 1 - 2_j Xi 



r(ai)---r(a„ + i) 1 



IaJv) 



1=1 



where A n = {x : xi > 0, Y%=i x i - !}• 

2 The pdf of the beta distribution with parameters (a, /?) is f(x;a,/3) = r(a)r(p) x °'~ 1 '0- ~ ^^"^(o.oo)^)- 
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An alternative approach following similar steps as for two and three-dimensions is developed 
in Appendix C. Note that (|2*5|) for a — > 2~ reduces to the multivariable Brownian motion with 
propagator 



m \ 2 "'l^iv-^ol 2 



lim P(\x N -x \,P;a) = h— e »0 . (29) 
The propagator (|25jl is the solution of the spatial fractional diffusion equation 



<9P(f,t) ft 



a-l 



_(_V^*P(£i) (30) 



9* (2m) 

with the condition P(x, i; x', t) = S(x — x') at equal times and P — > as \x — x'\ — »• oo. 



3 The Langrangian Approach 

Consider the action of the one- dimensional "free" particle in Minkowski space 

^ dt 



S 



M 



/(«) 









;(x) 2 ; 







(31) 



where the function /(a) has dimensions [E] 



[T]\ — 
[L] 



(with [E], [L], [T] representing the dimen- 



sions of energy, lenght and time respectively) and is given by 



(32) 



In f!32j) R a is a constant with dimenions [R a ] = [E] 1 "^ . The Lagrangian is invariant under time 
translations. Applying the transformation t — » t — tt together with a rescaling by (3 — tf — t± 
the action becomes 



Sm = /(a) U 



,1 



2( a -l) 



dt. 



(33) 



We decompose now the orbits x(t) of the point particle into a classical path x c i(t) and the 
deviations q(t) which are referred onwards as the quantum fluctuations of the particle orbit. 
The split and the boundary conditions x(t) satisfy are 

x(t) = x d {t) + q(t); x(0) = x{-l) = 0. (34) 

The classical path correponds to the solution of the Euler-Lagrange equation 

Xd (x d )^ = (35) 

with boundary conditions x c ;(0) = x in and x c [(— 1) = Xj. The solution of (|35j) with the 
prescribed conditions is 



(36) 
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The vanishing of the paths x(t),x c i(t) at the endpoints of the time interval [—1,0] forces the 
quantum fluctuations to obey the conditions q(0) = q{— 1) = 0. Inserting the decomposition 
into the action we have 



S 



M 



c-i r 



{id + q) Q_1 dt. 



(37) 



Using the binomial expansion (a + b) r = J2T=o ( ^ j ° r k b k with r, (k) be a positive real 



(nonnegative integer) number, 
we obtain 



r(r-l)-(r-k+l) 



and keeping only terms quadratic in g's 



*» = /(«> ft)"/ , 



a 



l 

at— 1 



a 



2(a-l) 2 



-,3 \ 



rft (38) 



The first term is dominated by the contribution along the classical path and the midterm 
vanishes identically after partial integration with respect to q(t) and using the equation of 
motion (|35j). The action then is written as 



5 



M 



/(«) 



[Xi - Xf 



+ 



a 



2-q 



2(a 



(qf 



dt. 



(39) 



The "free" particle amplitude in configuration space reads 



Sat 



(Xi,U)—>(xf,tf) 



e' 



(40) 



Expression (|4U|) implies a natural factorization of the amplitude into one factor steming from 
the contribution of the classical path and another arising from the quantum fluctuations. 

Let to = > ti > ■ ■ ■ > fx = — lbea partition of the time interval [—1,0] by points 
of subdivision t ,ti, ■ ■ ■ ,tjy. The quantum fluctuations are discretized according to qi = q(ti), 
i = 0, ■ ■ ■ , N with qo = qN = 0. The normalization factor of fluctuations occuring in (|4U|) 
will be evaluated by writing the functional integral as an infinite product of integrals over the 
discretized fields qu as follows 



Vq e 



lim 



(2Txh)~ N ( 2lxh 9i.Xi- x f ]a) 



f N-l 

n 

fc=i 



dq k e: 



r>Sjv 



(41) 



where g(xi-x f ; a) = f(a) (^Y 1 -^w {%i ~ Xf)^ 1 and SW = \g(xi-x f ;a)Y*=i (?* - qk-if ■ 
To diagonalize the discretized action S N we expand the quantum fluctuations in Fourier modes 



12 N ~ l fkrrm\ N ~ 1 
9* = v4E^ )r m =Y.OTr m] k = l,-,N-l. (42) 



N£i V N J 



m=l 



3 This expansion also appears in the continuum time limit by solving the Sturm-Liouvillc problem: q + Xq = 
0; q(Q) = q(— 1) = 0, with the unique solution q(t) = X)^Li r n -</2sm.(rnrt). 



8 



Substituting into the discretized action and employing the ortogonality relation J2m=i 0^0™ 
8j ; k we get a product of independent Gaussian integrals which can be calculated giving the result 



N-l 



ff f dq k e^-^U^f = ( - M -| 2 -L . (43) 
k=i J -°° \g(xi-x f ;a)J ViV 



Combining (j4T)j) . (|4*T|) and (}4*3*|) the amplitude is given by 

= — L=/(a)l \ . v/ " i . (x i -x / )^Ie^K^M^,(44) 

If we analytically continue /3 from Minkowski to Euclidean space using (3 = e - *^" -1 )^ the 
Euclidean amplitude becomes 

<x f ,tf\x i ,t i > = -/=/(«) 5 -^-(li-x^e riA ^ " .(45) 

Pe 

We will show now that this amplitude coincides with the long time limit of the propagator 
f!17|) at fixed position. We write it as 

POO 

P(\x N -x \,p;a) = — - e lz e- plzla dz (46) 

27l\XN — Xo\ J-oo 

where p = P . To get the best asymptotic approximation of (I4T1) with (3 — > oo we 

(2m) ^ \xpf—xo\ a 

apply the method of steepest descent. For this we write 

I(p; a) = / e iz - pza dz = [ e h{z) dz (47) 
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where 7 is the real axis in the complex 2-plane. We deform the original contour of integration 
onto 7' so that it takes a steepest descent path through the saddle point of h(z) avoiding 
singularities and branch cuts. Then on this contour, most of the contribution to the integral 
comes from the region around the saddle point, and possibly the endpoints. 

The location of the saddle points of h(z) is determined by setting the first derivative of h(z) 
to zero. Thus we find 

z ( a ) = l -^e^T^r). (48) 

(pa) 

Expanding h(z) in a Taylor series around z we have 

h(z) = h(z ) + ^(z-z ) 2 h"(z ) + --- (49) 
= h (z ) + ^(z-z ) 2 e i9 \h"(z Q )\ + --- (50) 

where e l6 is the phase of h"(zo). The linear part of the new path 7' is 

z - z = re<+ (51) 
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so that 



dz = e^dr. (52) 

Here e 1 ^ controls the direction of the path. We choose the direction of the path to decrease as 
rapidly as possible such that the function h(z) along this path has the form 

h(z) = h(z )-^\h H (z )\ + — (53) 
Expression (JH3|) coincides to (JHOj) provided that 

* - (54, 

Thus the integral I(p; a) becomes 

J(p; a) = e h ^ J°° e~^ h "^ 2 fe] dr. (55) 

Substituting in (jSHD M*o) = (£)^ (l - ±) e^(^r), ^"(^l = ( a _l) ( pa )^ , O = f (^f) 
and a = 2 ^+ 2 fc ^ , with fc G Z, we get 

P( arjv-a:o ,)9;a) = -=/(a) 3 — \ x " ~ x °l ^ e .(56) 

As a check one can prove that 

6 ~^~ (5?) 

since the binomial expansion becomes a perfect square and hence all higher order contributions 
vanish in this limit. 



4 Conclusion 

In the super-diffusion case 1 < a < 2, we have generalized the transition probability density of 
a particle governed by a symmetric and stable Levy flight in Euclidean n-dimensional space. 
Following the Hamiltonian approach the propagator is written in terms of the Fox function 
and its derivative in one- and three-dimensions correspondingly. In n — 2 and n > 3 the 
propagator is given in integral form using the Bessel functions. The Lagrangian approach on 
the other hand provides the asymptotic behaviour of the propagator in the long time limit at 
fixed position. The result produced by the Lagrangian method has also been crossed checked 
using the steepest descent applied to the exact result derived from the Hamiltonian approach. 
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Appendix A 



Definition 1 

Let X be a random vector on R n with probabiliy distribution /i and characteristic function 
uj. If fi m = //*-••* n denotes the m-fold convolution of [i with itself then we say that X is 
infinitely divisible if for each m G N there exist X ljm ,X 2i?Tl , • • ■ , X m m mutually independent 
and identically distributed random variables such that the sum X l m + X 2 , m + ■ ■ • + X m m has 
the same distribution as X. 

Hence if X i)m has distribution /i m and characteristic function oj m then /i = /z™ and ui = uj™. 
The Levy-Khintchine representation states that a probability measure fi on R n is infinitely 
divisible iff we can write the characteristc function uj{p) = E[e l<p,x> ] in the form exp(ip(p)) 
where 

$0n) = i <p,c> -- <p,Ap> + [ (e i<p ' x> - 1 - - <P \ X ^ ) M(dx) (A.l) 
2 j_rj y 1 + \ \x\\ 2 J 

where Rq := R n /{0}, c G R n , A is a symmetric nonnegative-definite n x n matrix (called the 
Gaussian covariance matrix) and M is a cr-finite Borel measure on Rq (called a Levy measure) 
such that 

mm{l, \ \x\\ 2 }M(dx) < oo. (A. 2) 

If A — then /i is said to be purely non-Gaussian. The triplet [c, A, M] is unique and will be 
called the generating triplet of the infinitely divisible distribution /i. 

A special case of infinitely divisible laws are the stable laws in n = 1. Suppose that 
X, Xi, • • • , X m denote mutually independent random variables with a common distribution F 
and S m = YT=x Xi- 
Definition 2 

The distribution F is stable if for each m G iV there exist constants c m > and 7 m G R such 
that 

S m = c m X + 7 m (A. 3) 

and F is not concetrated at one point. F is stable in the strict sense if 7 m = 0. 

The symbol = means that the distributions of S m and X differ by location and scale pa- 
rameters. The norming constants are of the form c m = with 1 < a < 2 and the constant 
a is called characteristic exponent of F or index of the stable law. 

If we consider A = and the Levy measure to be of the form M(r, oo) = pCr~ a and 
M(— oo, — r) = qCr~ a where r, C > 0, p, q > 0, p + q = 1 for some < a < 2 then fi has 
characteristic function given by 

1 -ibsgn(p)tim(?f)\ ; a^l 
^(p) = «7P-s m , r, , 2 / \ i /i .J , (A.4) 




1 + -sgn(p) ]n(\p\) 



a = 1 



where b = p — q, c = C r( f~^ cos( 2 ^) with — 1 < 6 < 1, c > and j G R. The parameters 
a, 6, c, 7 characterize the asymptotic behaviour, the skewness, the scale and the location of the 
peak of the stable distribution respectively. Moreover the collection (a,b, c, 7) is called the 
stable law parameters. A stable law generated by (0,6,0,7) is often denoted by 5*0(6,0,7). In 
the present work we set the degree of asymmetry to zero (6 = 0). 



11 



The ath order Weyl derivatives on the infinite axis are defined as [H] 



def 



d r 



dx m 
( d X 



rl-{a\ 

+ 

d< 



d 



Tim — a) \ dx 



fit) 



(x - t) a ~ m+1 



dt, 



dx r 



Tim — a) \dx 



" {l-~ {a} f(x] 

fit) 
(t - x) a - m+1 ' 



■dt. 



(A.5) 



where a = [a] + {a} with [a], {a} standing for the integral and fractional part (0 < {a} < 1) of 
the real number a. Also m = [a ] + l and V%f{I± w ) are the left- and right-handed fractional 
derivatives (integrals) and V the Euler's gamma function. It can be shown that for / e C^°(Q), 
Q C R the Fourier transforms of (|A.5|) satisfy for < a < 2 

T (Vlf{x)) = {TivT T {f{x)) = (TipT fiv) (A.6) 

where !F{f{x)) = f{p) = J R e tpx f(x)dx and (Tw) a — \p\ a e ¥ ' l ^ as9n(j '^ . One can prove that 



V*Vlf{x) =T~ 



-ip) Q (ip) Q f(p) 



WKp) 



2\ 2 



fix) 



(A.7) 



where V 2 is the one dimensional Laplacian. 

Proposition 1 

Let fix) = e tp ' x where p, x e Then 



(A.8) 



(-V^/(x) = (p 2 )*/(x) 

where p 2 = Y%=i Pi an d Pi £ -R- 
Proo/ 

Define the translation operator of a function /(x) by 

(T n f)(x) = f(x-h) (A.9) 
and the finite difference of order I of f{x) with step /i and center at the point x, Ai, by 



(Al/)(f) = (/-T K )V(xl = E( / A; ) 



(-!)*/(*- fc/o 



(A.10) 



where / is the identity operator. It can be shown that the Fourier transform of the hypersingular 
integral 



D a f 



R n \y 



n+a 



d n y 



is given by 



F{D a f{x)) = \p\ a d n i{a)F{f{x)) 



where the constant d n i(a) is 



d n Aot) 



d n u 



'1-e 



ip-u\l 



III 



n+a 



(A.H) 



(A.12) 



(A.13) 
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In (IA.13|) p is the unit vector in the direction of p. Using as f(x) = e l $' x we have 

1 



■V 



2\ 2 ip-x 



d n ,l( a ) 
1 



a ip-x 



D a e 



d n ,i(a) Js* \y\ n + a ^ Q \ 
1 ** f ^d-e^7 



\y\ 



n+a 



(A.14) 



Appendix B 

The class of Fox or H-functions [2j comprises a large class of special functions known in math- 
ematical physics. It is defined in terms of the Mellin-Barnes path integral 



rjm.ni ' J\ Tjm,n 



(ai, Ai), (a 2 , A 2 ), ■ ■ ■ , (a p , Ap) 
(6i,Si),(6 2jJ B 2 ),... J (6„S,) 



^-J L X(sy d s (B.l) 



with the integral density 



n? 1= m+l ^(l — fajj + Bj l s) rii 2 =n+l r(a« 2 — Aj 2 S, 



The H-function possesses a number of interesting properties jH] which we list the 
employed in our presentation. 

• If k > then, 



(B.2) 



ones we 



2C B (*) = 



2 



(ai, Ai), (d 2 , A 2 ), 
(b 1 ,B 1 ),(b 2 ,B 2 ),---,(b q ,B 



(B.3) 



Under Fourier cosine transformation, the H-function transforms as 



Tjm,n 

P.9 



(dp, A p ) 
(6,,^) 



cos(foc)da; = 



(l-6 g ,B 9 ), (1,|) 

(i,i),(i-dp,Ap),(i,|; 



• (B.4) 



The function z b e 2 can be represented in terms of the H-function as 



1,0 



rr±,u 
-"0,1 



(0,0) 
(6,1) 



(B.5) 
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Appendix C 

Using spherical coordinates in (|22|) one has 



P(\x N -x \,P;a) = j o dpNp^e-x^Ftp^lxN-Zoln) (C.l) 



where 



tn-2 



F(p N ,\x N -x \,n)= l^ilyll J* d6 k sin*- 1 e^^dO,^ sm n ~ 2 e^e^"-^™ 9 ^.^) 
Performing the integrations over the angles the integral representation of the propagator is 
P(\x N - SolPta) = (^) 2 ^.^^-i /"dP* ph-^ N) J^i{^-\x N - x \). (C.3) 
The calculation of the integrals over the angles is based on the formulas [0] 

Jo Vr(±(m + 2))' 

q d6 sin"W ccose = | 1 ^x(l-x 2 ) :2 f i e icx = 2 1 -f7r^|c-f J f (c) 

hl-aPY-icos^dx = ^ f-Vr(z/+i)J,(a) 
Jo 2 \a/ 2 

/■°° „ _ aa: 2 r(z/ + /i + 12 7? 2 

Jo j3a-2Y[y + 1) 22 4a 

2 2:E - 1 r(x)r(x + = v^Fr(2x). (c.4) 

In (IC4J) is the Bessel function of the first kind of order v defined by the equations [TO] 

J ^ - id) 

z \ v 00 z 2 ^ 

2) ^ ( " 1)k »HHl) (C,5) 



or in a more suitable form by 



where n is a natural number. For completeness we give the expicit expressions for J and Ji 

2 

used for the propagators in two and three dimensions. They are 



„2k 

M*) = E(-!) 

fc=0 



2 2k (k\) 2 



Ji(z) = J—smz. (C.7) 
V nz 



2 
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